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Processes  exist  today  which  will  produce  oxide  layers  on  silicon  semiconduc- 
tor devices  with  total gonizing  surface  effects  sufficiently  suppressed  to 
satisfy  military  needs.  There  remains  the  problem  of  assuring  that  the 
devices  produced  by  these  processes  on  a production  line  are  indeed  hard  to 
surface  effects,  i.e.,  an  unintended  variation  in  processing  has  not 
reduced  hardness.  Some  sort  of  acceptance  test,  called  a hardness  assurance 
screen,  is  needed.  To  identify  possibilities  for  these  tests,  the  litera- 
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“ ture  was  searched,  analyses  were  performed,  and  special  wooden  ball  molecular 
models  were  studied.  The  current  model  of  oxide  layer  growth  and  the  rela- 
tions of  such  growth  to  total  dose  susceptibility  were  gleaned  from  the 
literature.  This  model  was  extended  by  the  study  as  follows:  Qss'  i s caused 

by  an  oxygen  atom  which  is  bonded  to  thr^e  interface  silicon  atoms.  This 
produces  an  interface  donor  state  whose  energy  is  above  the  silicon  conduc- 
tion bano.  N’$  is  primarily  caused  by  a Si -0-Si  bond  in  which  the  oxygen 

has  been  removed  leaving  a stretched  Si -Si  bond.  This  produces  two  interface 
states;  a donor  level  near  the  valence  band  and  an  acceptor  level  near  the 
conduction  band..*  N„  near  the  band  edges  has  the  same  cause  as  N . The 

oxygen  is  removegVrom  the  Si -0-Si  bond  by  direct  action  of  the  radiation 
and  by  the  attractiye  action  of  the  positively  charged  ions  which  produce  Q 
Qr  is  a result  of  the.  holes  which  are  produced  by  the  ionizing  radiation 

being  trapped  on  silicon  excess  centers.  Silicon  excess  centers  are  centers 
where  one  oxygen  is  bonded  to  three  silicon  atoms  in  the  bulk  of  the  SiOj. 

The  excess  centers  are  distributed  according  to  the  form 

A„ . = K(e"ax  -'ea(x-2t))/(l  + e'2at) 

\ 

The  susceptibility  to  Qr  is  proportional  to  Qss  measured  after  oxidation, 

before  annealing.  Internal  electric  fields  are  important  in  oxide  growth 
processes  and  the  growth  can  be  influenced  by  anything  which  alters  the 
field.  Finally,  the  extended  model  was  examined  for  suggested  hardness 
assurance  screens  and  suggested  experiments  to  confirm  the  extentions. 
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SECTION  I 
INTRODUCTION 


The  problem  addressed  by  this  report  stems  from  efforts  to  harden  silicon 
semiconductor  devices  to  surface  or  total  dose  effects.  Two  surface  effects 
are  produced  by  the  action  of  ionizing  radiation  on  the  silicon  dioxide  layers 
common  to  silicon  semiconductor  devices.  One  of  these  effects  it  the  apparent 
accumulation  of  positive  space  charge  (Qr).  The  second  effect  is  the  produc- 
tion of  electrically  active  defect  centers  located  at  the  interface  between 
the  silicon  and  the  silicon  dioxide  layer.  These  centers  are  called  interface 
states  (Nr).  During  the  last  few  years,  empirical  efforts  have  been  made  to 
minimize  these  effects  by  varying  the  growth  parameters  (temperature,  thick- 
ness, oxidizing  gas,  etc.)  of  the  silicon  dioxide  layer  or  by  adding  impur- 
ities (aluminum,  chromium)  to  the  layer.  Processes  exist  today  which  will 
produce  oxide  layers  on  silicon  semiconductor  devices  with  these  surface 
effects  sufficiently  suppressed  to  satisfy  military  needs.  There  remains  the 
problem  of  assuring  that  the  devices  produced  by  these  processes  on  a produc- 
tion line  are  indeed  hard  to  surface  effects,  i.e.,  an  unintentional  variation 
in  processing  has  not  reduced  hardness.  Thus,  some  sort  of  acceptance  test, 
called  a hardness  assurance  screen,  is  needed. 

Work  in  the  past  has  shown  that  surface  effects  could  be  annealed  out  of 
the  oxide  layer  by  baking  (e.g.,  300°C  for  24  hours).  One  successful  screen- 
ing procedure  consists  of  first  irradiating  a device,  thus  determining  its 
radiation  susceptibility,  and  then  annealing  the  induced  effects  (reference  1). 
This  procedure  is  very  Involved  and  costly.  A simple  electrical  measurement 
involving  no  special  ambient  would  be  much  less  costly. 

Two  approaches  are  possible.  One  approach  is  to  experimentally  measure 
all  measurable  properties,  determine  the  radiation  susceptibility  by  a radia- 
tion test,  then  look  for  correlations.  The  second  approach  (which  is  the 
subject  of  this  report)  is  to  examine  the  most  detailed  model  of  the  resoon- 
sible  basic  mechanisms  to  see  If  any  correlations  between  measureables  and 
hardness  are  suggested. 

The  procedure  of  examining  the  most  detailed  model  available  normally  will 
yield  no  new  insights  unless  the  model  is  first  extended  to  include  some  new 
ideas.  Then,  such  a model  may  suggest  a screen,  and  it  may  suggest  variations 
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in  processing  which  will  produce  harder  device  oxide  layers.  Here,  we  extended 
the  model  by  Including  the  electric  field  present  during  growth,  so  that  the 
currently  accepted  models  of  thermal  oxide  growth  were  expanded  to  include 
anodic  oxide  growth.  Also,  molecular  models  of  the  silicon  lattice  were  used 
to  study  the  oxidation  process  atom  by  atom  and  to  explore  the  structure  of 
interface  states. 

There  are  two  surface  properties  which  can  be  measured  by  simple  electri 
cal  methods.  One,  the  surface  state  charge,  Qss  (normally  given  in  electronic 
charges/cm2),  produces  the  same  effects  in  the  silicon  as  a sheet  of  electro- 
static charge  and  seems  unaffected  by  any  surface  potentials.  Two,  the  inter- 
face state  density,  Ns$  (given  in  states/cm2/ eV ) , is  a spectrum  of  energy 
levels  located  in  the  forbidden  gap  of  silicon  which  can  trap  and/or  recom- 
bine electron  and  hole  carriers.  The  post-irradiation  properties  Qr  and 
are  the  radiation  induced  properties  detected  by  the  same  methods  as  Qss  and 

Nss- 

In  the  following  sections,  the  insight  obtained  from  extending  the  model 
are  used  to  develop  hypotheses  about  Qss  and  N$s.  We  propose  that  the  rate 
that  silicon  Is  stripped  off  the  silicon  lattice  during  oxide  growth  is  pro- 
portional to  a function  of  temperature  and  oxygen  partial  pressure  times  Qs$. 
For  increases  in  rate  due  to  Increases  in  temperature  or  electric  field,  Q$s 
decreases.  For  increases  in  rate  due  to  increases  in  oxidation  nucleation 
sites,  Q$s  increases.  It  Is  hypothesized  that  Nss  is  a direct  measure  of  the 
oxygen  chemical  activity  at  the  silicon-silicon  dioxide  interface.  We  relate 
the  radiation  susceptibility  of  the  oxide  layer  to  the  amount  of  excess  sili- 
con In  the  oxide.  The  amount  of  excess  silicon  is  determined  by  the  balance 
between  the  rate  of  silicon  introduction  and  the  rate  of  oxidation  of  the 
extra  silicon,  i.e.,  the  oxygen  activity.  Thus,  we  propose  that  a function  of 
temperature,  Qs$  and  Nss>  may  prove  to  be  the  desired  hardness  assurance 
screen  and  that  variations  of  either  Q$s  and  N$s  are  cause  for  concern. 
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SECTION  II 

SUMMARY  OF  THE  PRESENT  MODEL 

In  March  1967,  Deal,  et  al.  (reference  2)  suggested  the  following  model 
for  thermal  oxide  growth: 

a.  At  the  oxygen-silicon  dioxide  interface,  oxygen  enters  the  SIC2  in_ 
creasing  the  stoichiometric  proportion  of  oxygen  in  the  SI O2.  above  the  normal 
two  oxygen  atoms  for  each  silicon  atom.  The  excess  oxygen  cannot  be  uni- 
formly distributed  on  an  atomic  scale  but  must  be  located  in  discrete  centers. 
The  structure  of  the  centers  can  be  an  interstitial  oxygen  or  a single  dang- 
ling oxygen  bond,  called  a non-bridging  oxygen. 

b.  At  the  silicon-silicon  dioxide  interface,  silicon  interacts  with  sili- 
con dioxide  to  form  silicon  excess  centers.  A silicon  excess  center  car.  be  an 
oxygen  vacancy,  a single  dangling  silicon  bond,  or  an  oxygen  atom  bonded  to 
three  silicon  atoms  instead  of  the  normal  two. 

c.  The  above  centers  diffuse  into  the  bulk  Si02.  Since  the  excess  oxygen 
center  diffuses  more  rapidly  than  the  excess  silicon  center,  the  two  recombine 

O 

near  the  silicon  interface,  i.e.,  within  100  to  200  A. 

d.  After  growth,  the  oxide  contains  oxygen  excess  centers  except  for  a 

O 

region  100  to  200  A from  the  silicon  interface  which  contains  silicon  excess 
centers. 

This  model  could  explain  the  production  of  positive  space  charge  (Q  ).  It 

o r 

is  known  that  this  charge  is  trapped  within  200  A of  the  silicon  interface 
(reference  2).  A commonly  proposed  model  is  that  this  positive  space  charge 
is  trapped  In  the  excess  silicon  centers,  item  b.  of  Deal's  model.  There  is 
no  generallynaccepted  model  of  what  causes  the  increase  in  interface  states 

(Nf). 
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SECTION  III 

EXTENSIONS  OF  THE  MODEL 


In  this  section  the  model  is  extended  to  explain  additional  data  about 
Qss>  The  electro-chemical  nature  of  the  oxide  film,  and  the  lattice  structure 
of  the  silicon  at  the  interface. 

1.  Surface  Charge  Qss 

Deal  (reference  2)  proposed  that  (1  was  a volume  charge  distributed  in 
the  100  to  200  A of  oxide  next  to  the  silicon  interface.  However,  Powell 
(reference  3)  performed  an  experiment  which  is  best  explained  if  (3  is  with- 

o SS 

in  20  A of  the  interface.  If  this  is  true,  Qss  cannot  be  a volume  charge. 
Powell's  experiment  suggests  that  Qss  is  an  interface  state  and  is  not  lo- 
cated on  the  same  centers  as  Q . 

r 

If  Qss  were  a volume  charge,  it  would  be  trapped  on  the  same  centers  as 
the  radiation  charge  Qr<  Thus,  the  size  of  Qss  would  depend  directly  on  the 
number  of  centers  as  would  the  size  of  Qr<  There  would  then  be  a direct  cor- 
relation between  Q$s  and  Qr.  A correlation  as  direct  as  this  has  not  been 
generally  observed. 

If  Qss  were  a volume  charge  and  distributed  on  the  excess  centers,  then 
because  positively  charged  excess  centers  are  unfilled  electron  states. 


«,sa  • 4X 


(1) 


where  aX  Is  the  thickness  of  the  layer  (100  A) 

Nes  Is  the  average  density  of  excess  centers 
f is  the  Fermi  factor 
(1-f)  is  approximately  exp  [(E.  - E^.)/kT j 
E^  Is  the  energy  level  of  the  excess  center 
Ef  is  the  Fermi  level  of  electrons 

In  the  silicon,  at  oxide  growing  temperatures  (1000°C),  the  Fermi  level  E- 
must  be  in  the  center  of  the  gap  as  shown  in  figure  1.  In  going  from  the 
silicon  to  the  silicon  dioxide,  the  electrostatic  potential  must  change  by 
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SILICA  CONDUCTION  BANO 


OXYGEN  EXCESS  CENTER  ACCEPTOR  LEYEL 


1.1  8V 


Figure  1.  Band  Diagram  of  S i - S i 0 2 Interface. 


the  amount  of  the  electrochemical  potential,  1.616  V at  1000°C  (table  1). 

The  value  of  is  about  1 eV  above  the  Si02  valence  band  (reference  4)  which 
is  in  turn  3.8  eV  below  the  valence  band  of  silicon  (reference  3).  Then  (1-f) 
= 8 x 10"9  and  using  Qss  = 1011  charges/cm2,  Ngs  = 1.2  x 1024,  which  is  about 
102  times  the  atomic  density  and  is  quite  impossible.  Temperature  correc- 
tions to  take  the  above  room  temperature  energy  levels  to  luuO°C  are  about  5 
percent  (reference  5),^hich  would  not  change  the  result. 

Thus,  the  hypothesis  that  (T  is  a volume  charge  and  is  distributed  in 

o S S 

the  100  to  200  A band  near  the  silicon  interface  as  proposed  by  Deal  does  not 
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Table  1 . 

ELECTROMECHANICAL  POTENTIALS  FOR  SILICON  AND  OXYGEN  REACTIONS 

AFt(volts) 


Temp 

°C 

2Si  + 02 
-►2SiO* 

Si  + 02 
-Si02** 

Si  + Si02 
+SiO 

2Si0  + 02  u 
+S102 

Si0“ESi02)/kT 

25 

-2.043 

-2.058 

+.0153 

-2.074 

0.55 

65 

2.031 

2.040 

+.0091 

-2.049 

0.73 

300 

1.956 

1.932 

-.0246 

1.917 

1 .64 

400 

1.924 

1.886 

-.0376 

1.848 

1.91 

500 

1.890 

1.840 

-.0497 

1.791 

2.11 

600 

1.846 

1.795 

-.0612 

1.734 

2.25 

800 

1.787 

1.705 

-.0819 

1.623 

2.42 

900 

1.752 

1 .661 

-.0912 

1.569 

2.46 

1000 

1.716 

1.616 

-.1000 

1.516 

2.49 

1100 

1.686 

1.572 

-.1081 

1.464 

2.49 

1200 

1.644 

1 .529 

-.1157 

1.413 

2.49 

*AFt 

= -97709.  +3.84  TlogT 

+ ( . 1 6E-3 ) T2  - 

C.295E  + SJT’1 

+ w.45T,  Bull e- 

tin 
of  - 

542,  U.S.  Bureau  of  Ml 
-.76,619  cal  at  1880°C. 

nes,  1954,  corrected  using  heat  of  sublimation 

**AFt 

= -203,380. 

- .43TlnT 

- (1.06E-3)T2  + 

( . 98E  + sjr1 

+ 47.01T. 

fit  the 

data  and  can 

be  rejected  in  favor  of 

the  hypothesis 

that  Qss  is  an 

interface  state.  This  extensi 
with  an  energy  level  above  the 

2.  Constraints  Imposed  by 

on  of  the  model  says  that  Qss 
conduction  band. 

the  Lattice  Structure 

is  a donor  state 

The  silicon  crystal  lattice  is  by  definition  a regular  structure  with 
short  and  long  range  order  and  with  certain  separations  between  atoms.  Ob- 
viously, this  structure  must  impose  limitations  on  the  structure  of  the  inter- 
face between  the  silicon  lattice  and  the  amorphous  silicon  dioxide  surface 
layer.  The  purpose  of  this  section  is  to  theoretically  explore,  using  wooden 
ball  molecular  models,  possible  surface  structures  and  the  attendant  implica- 
tions. The  study  is  restricted  to  considering  only  structures  involving  sili- 
con and  oxygen  atoms.  With  this  restriction,  only  three  bonding  configura- 
tions are  probable.  First,  the  oxygen  can  have  one  bond  with  the  silicon 
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lattice  and  either  one  or  none  with  the  amorphous  Si02.  Second,  the  oxygen 
can  have  two  bonds  with  the  silicon  lattice.  The  third  configuration  con- 
sists of  triply  coordinated  oxygen  with  three  bonds  to  three  different  silicon 
atoms . 

In  the  following  discussion,  two  bonding  structures  will  be  defined. 
Standard  quantum  mechanical  considerations  of  the  way  the  energy  levels  of 
these  structures  must  split  will  qualify  them  as  candidates  for  the  structure 
of  the  defects  which  cause  Qss>  the  fixed  interface  charge  resulting  from 
oxidation,  and  Ns$,  the  interface  state  density  following  oxidation.  To  find 
experimental  confirmation  that  these  structures  are  indeed  Qss  and  Nss>  the 
Implications  of  the  hypothesis  are  explored.  By  studying  the  way  they  are 
created,  Nss  is  related  to  oxygen  activity  during  oxide  growth  and  Qs$  to  the 
rate  of  oxidation.  By  assuming  that  stoichiometrically  excess  silicon  intro- 
duced during  growth  is  the  cause  of  the  hole  traps  which  produce  flat  band 
shift,  a complex  formula  predicting  effective  hole  trap  density  as  a function 
of  Qss»  Nss>  temperature,  and  oxide  thickness  is  derived.  The  formula  is  not 
inconsistent  with  known  experimental  data,  and  it  suggests  new  experiments. 

Whe.  cxyge.  forms  a single  bond  with  the  silicon  lattice,  it  would  nor- 
mally have  a second  bond  with  the  amorphous  Si02.  However,  attempts  to  attach 
Si02  tetrahedra  to  the  silicon  lattice  yielded  a negative  result.  Basically, 
the  lattice  constraints  of  the  crystal  forms  of  silica  are  so  different  from 
silicon  that  order  cannot  be  transferred  across  the  interface.  This  situa- 
tion is  shown  in  figure  2,  which  actually,  is  very  unreal.  The  <1 11 > surface 
is  shown  bonded  to  silica  tetrahedra.  In  silica  the  Si -0-Si  bond  is  from  3 to 

O 

3.1  A long.  The  <1 1 1 > surface  (see  figure  3)  has  silicon  atoms,  with  one 

o 

dangling  bond  per  atom,  in  equilateral  triangles  3.84  A on  side,  i.e.,  in 
figure  1,  the  silica  tetrahedra  wouldn't  be  tetrahedra.  If  order  is  not 
transferred  across  the  interface,  then  there  must  be  frequent  defects.  The 
knowledge  that  the  oxide  reduces  1015  dangling  bonds  of  the  cleaved  lattice  to 
between  1010  and  10u  interface  states  (reference  2)  indicates  that  such 
single  bonds  are  uncommon. 

The  possibility  of  oxygen  making  a single  bond  (nonbridging  oxygen)  with 
the  surface  cannot  be  discarded  on  theoretical  grounds.  However,  such  a bond 
would  produce  an  interface  accepter  level  near  or  below  the  valence  band. 

These  have  not  been  observed. 
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Figure  2.  Unsuccessful  Attempt  to  Place  Silica  Tetrahedra  on  the  Silicon  <111>  Plane 
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Considering  the  single  bond  cases  results  in  the  conclusion  that,  bridging 
bonds  between  the  lattice  and  the  oxide  cannot  exist  in  any  great  number 
(i.e.,  > 1010  - 1 0 1 1/cm2 , nonbridging  oxygen  would  produce  an  unobserved 
energy  level). 

Next,  consider  the  double  lattice  to  oxygen  bond  case.  Here  oxygen  forms 
bonds  with  two  different  Si  atoms  in  the  lattice.  In  figure  3,  the  available 
Si  bonds  are  shown  for  the  <100>,  <11 0>  and  <111 > planes.  Surprising  is  the 
fact  that  the  broken  bond  density,  a,  is  the  largest  for  <1 OC > . This  is  the 
reverse  of  the  occurrence  of  Qss,  the  surface  charge  density,  and  N$s>  the 
interface  state  density.  Figure  4 shows  the  simple  way  that  oxygen  atoms  can 
be  double  bonded  to  the  silicon  lattice  in  the  <100>  case.  Because  of  the  way 
the  wooden  balls  were  drilled,  figure  4 shows  a Si -0-Si  bond  angle.  Silicon 

O 

atoms  in  the  S i -0-Si  bonds  of  figure  3 are  3.84  A apart.  This  bond  is  3 to 

O 

3.1  A in  silica.  Accordingly,  this  bond  is  probably  straight,  as  it  is  in 
B - quartz. 

In  the  <11 0>  and  <11 1>  cases,  the  Si -0-Si  bonds  can  be  made  in  more  than 

one  way.  The  ways  are  illustrated  in  figure  5 <11 0>  and  figure  6 <11 1 > . As 

can  be  seen  in  figure  3,  the  atoms  of  the  <11 0>  surface  are  in  chains  as 

0 

shown  in  figure  5.  These  chains  are  separated  by  5.43  A centerline  to  center- 

O 

line.  The  atoms  in  different  chains  are  at  least  4.5  A apart.  Thus,  the 
passivation*  of  the  chains  can  proceed  independently  of  each  other.  In  figure 
5,  locations  for  the  passivating  oxygen  location  atoms  are  numbered.  To  passi- 
vate the  chain,  only  every  other  oxygen  location  on  each  side  of  the  chain  is 
filled;  i.e.,  all  the  1's  or  all  the  4's  and  all  the  3 ' s or  all  the  2's. 
Clearly,  this  is  a more  complicated  passivation  structure  and  defects  are  more 
probable  than  in  the  <100>  case. 

From  figure  3,  it  can  be  seen  that  the  surface  atoms  on  the  <111>  surface 
have  an  equilateral  triangle  pattern  as  represented  by  the  small  circles  in 
figure  6.  In  figure  6,  the  possible  oxygen  sites  are  large  circles.  Only  a 
sixth  of  the  circles  actually  have  oxygen  atoms  in  th^m  when  the  surface  is 
fully  passivated.  In  figure  6,  six  oxygen  sites  about  one  atom  have  been 
numbered.  Other  sites  have  numbers  with  a line  over  them,  e.g.,  3".  The 

*The  S1-0-S1  bond  has  the  ionization  energies  of  silica,  i.e.,  very  far  above 
and  below  the  band  gap  of  silicon.  Because  the  electrons  in  this  bond  are 
accordingly  made  inactive,  the  making  of  these  bonds  is  called  passivation. 
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Broken  Bond  Patterns  of  the  Si  Interface 
a is  the  Broken  Bond  Density. 


Figure  3 
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Figure  6.  Oxygen  Sites  on  the  <11 1>  Surface. 

overlined  numbers  mean  that  if  there  is  an  oxygen  atom  in  the  numbered  site, 
there  cannot  be  one  in  this  site  and  vice  versa.  Clearly,  this  <111>  face  has 
a still  more  complicated  passivation  structure  than  the  <1 1 0>  face  and  defects 
would  be  still  more  probable. 

The  above  discussion  shows  the  way  tha't  the  dangling,  energetically  un- 
favorable, broken  bonds  of  a perfectly  cleaved  silicon  surface  can  be  passi- 
vated with  stable  Si-O-Si  double  bonds.  Because  no  surface  is  perfectly 
cleaved  and  because  entropy  forbids  a perfectly  orderly  passivation  of  such 
surfaces,  there  will  be  defects  in  this  order.  Such  a defect  is  a pair  of 
dangling  bonds  with  no  oxygen  atom.  Specifically,  one  of  the  oxygens  in  fig- 
ure 4 could  be  removed  leaving  two  unpassivated  bonds.  A similar  situation 
occurs  in  bulk  silicon.  In  figure  7 the  S i A center  (reference  7)  is  shown. 

The  Si A center  is  a defect  in  the  bulk  caused  by  a substitutional  oxygen.  The 
levels  in  this  center  are  produced  by  the  stretched  Si-Si  bond.  Because  the 
silicon  atoms  are  further  apart,  the  splitting  between  the  bonding  (valence 
band)  and  the  anti-bonding  (conduction  band)  is  less.  Thus,  two  levels  are 
created  In  the  forbidden  gap.  Similar  levels  would  be  produced  by  such  a 
stretched  bond  on  the  surface.  The  above  yields  a model  for  observed  states 
near  the  band  edges  (references  4,  5,  and  6).  Also,  as  shown  in  figure  7 (v^) 
the  energy  level  of  a single  dangling  bond  which  is  now  unsplit,  would  be  in 
the  center  of  the  band  gap.  For  the  remainder  of  this  report,  let  Nss  denote 
the  levels  near  the  band  edges  only. 
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Figure  7.  The  S1A  Center  (reference  7). 

The  number  of  these  states  (Nss)  which  are  not  passivated  by  an  oxygen  atom 
will  be  determined  as  thermodynamically  required  by  the  law  of  mass  action,  i.e 

Nss  A0  = F(T)e-Es/kT 


(2) 
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where  Aq  is  the  oxygen  activity  (*  oxygen  atom  concentration  at  the  silicon 
interface) 

T is  absolute  temperature 

F(T)  Is  some  varying  function  of  temperature  and  has  a radically  differ- 
ent value  for  different  Si  lattice  planes 

E is  the  activation  energy  for  separating  the  oxygen  from  the  pair 
s 

of  dangling  bonds 

Accordingly,  a lack  of  oxygen  near  the  interface  produces  N$$,  and  N$s  is  in- 
versely proportional  to  oxygen  partial  pressure. 

Conditions  which  permit  triply  bonded  oxygen  can  be  found  by  exploring 
the  way  the  silicon  surface  oxidizes.  Each  plane  must  be  penetrated  by  re- 
moving a silicon  atom.  Then  succeeding  atoms  in  the  plane  are  removed,  spread- 
ing in  circular  waves,  peeling  the  plane  away  like  the  leaves  of  a magazine  in 
a fire.  The  nature  of  this  exposed  plane  edge  is  explored  in  figure  8.  In 
the  <100>  case  in  figure  8,  notice  the  three  exposed  bonds.  In  the  oxidizing 
process  the  oxygen  atoms  take  the  position  of  the  silicon  atoms  to  which  they 
are  bonded.  The  displaced  silicon  atoms  enter  the  oxide  and  the  starting 
structure  is  repeated.  Thus,  the  plane  is  progressively  stripped  off.  In  the 
<110>  case,  planes  are  stripped  off  by  nipping  off  the  end  of  the  chain  (fig- 
ure 6).  The  process  of  stripping  in  each  case  produces  a local  region  much 
more  oxygen  deficient  than  other  areas  of  the  interface.  Imperfect  passiva- 
tion is  likely  and  pairs  of  unpassivated  bonds  do  not  have  time  to  diffuse  to 
adjacent  sites.  Centers  of  the  type  shown  in  figure  9 are  likely  to  occur. 
Figure  8 shows  how  this  center  occurs  in  each  of  the  cases.  The  center  has 
two  silicon  atoms  with  their  broken  bonds  well  screened  by  the  oxygen.  The 
third  silicon  has  its  fourth  bond  partially  screened  by  the  same  oxygen.  This 
situation,  studied  in  Si02  glass  by  Weyl  (reference  8)  results  in  a center 
with  a very  loosely  bound  electron  where  a net  positive  charge  is  energeti- 
cally favorable.  At  the  interface  such  a structure  with  a positive  charge  and 
the  neutralizing  electron  orbiting  in  the  silicon  conduction  band  could  be  the 
source  of  the  positive  fixed  charge  Qss-  The  resulting  level  would  be  above 
the  conduction  band  edge  where  the  Fermi  level  would  never  reach  it,  thus  a 
fixed  charge. 
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Next,  consider  the  way  successive  planes  of  silicon  are  removed.  Given 
the  assumption  that  the  silicon  oxidizes  more  readily  at  the  plane  edges  than 
on  the  surface  of  the  plane,  the  condition  discussed  above  results  in  a model 
of  oxide  growth.  The  oxidation  of  each  successive  plane  must  begin  with  the 
puncture  of  the  plane  surface.  The  plane  is  then  quickly  oxidized  off,  un- 
covering the  next  plane.  A puncture  can  occur  just  as  readily  on  a plane 
which  is  already  punctured  as  one  which  is  not.  If  this  puncturing  is  a 
random  process,  then  the  probability  of  a puncture  must  be  proportional  to 
the  exposed  area.  This  situation  can  be  visualized  as  rain  falling  in  a pool 
of  water;  the  rain  drops  are  the  punctures  and  the  expanding  wavelets  are  the 
oxidizing  edges.  The  analogy  is  not  perfect,  in  that  in  the  pool  of  water, 
two  wavelets  flowing  together  do  not  disappear.  However,  it  illustrates  that 
there  are  two  important  independent  parameters,  the  rate  of  punctures  and  the 
wave  velocity  of  the  oxidizing  edge.  Each  must  obey  the  Arrhenius  form, 

R = Nve"Ea/kT  (3) 

where  N Is  the  number  of  available  sites 

v is  the  vibration  frequency  OO13  sec-1 
kT  is  the  thermal  energy 

E.  is  an  activation  energy  which  can  be  modified  by  electric  fields 

a 

The  structure  shown  in  figures  8 and  9 requires  that  Q$s  be  proportional 
to  the  total  of  all  oxidizing  edge  arc  lengths  remaining  after  cooling.  Then 
Qss  increases  with  puncture  rate  and  decreases  with  oxidizing  wave  velocity. 

It  is  known  that  Q$s  decreases  with  temperature  (reference  2).  This  effect 
can  be  produced  if  the  activation  energy  of  the  wave  velocity  is  greater  than 
the  activation  energy  of  the  puncture  rate.  Observe  that  each  puncture  (which 
oxidizes  one  silicon  atom)  sets  off  a wave  which  oxidizes  many  silicon  atoms. 
Thus,  most  of  the  silicon  atoms  come  from  the  wave  front.  Since  the  total 
wave  front  arc  length  is  proportional  to  Qss , then  from  equation  (2),  the  rate 
that  silicon  atoms  enter  the  oxide  is 

R ■ R0  Qsse*P  <-EssAT)  (4) 
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where  Rq  is  a function  of  temperature  equal  to  Nv/Qss 

Ess  is  the  activation  energy  for  the  wave  front. 

There  is  no  reason  to  believe  that  the  ratio  of  N to  Q$s  does  not  vary  with 
temperature  so  RQ  contains  temperature  components  from  N/Qss  and  v.  RQ  will 
also  contain  Aq  as  a factor,  as  shown  in  the  next  section  (3). 

Thus  the  lattice  produces  at  the  edge  of  a plane  a region  where  three 
bonds  must  be  made  with  oxygen  rather  than  the  energetically  more  favorable 
two.  This  produces  a model  relating  Q$s  to  growth  rate. 

3.  The  Electric  Field  at  the  Interface 

In  1962  Jorgensen  (reference  9)  showed  the  electrolytic  nature  of  the 
Si02/Si  system  by  making  a reversible  Pt/Si02/Si  electrolytic  cell.  Hot,  very 
pure  Si 02  must  then  be  an  electrolyte,  i.e.,  be  ionized.  Here,  the  possible 
nature  of  the  ions  in  Si02  will  be  explored,  as  well  as  the  effects  of  the 
electrolytic  nature  has  on  the  growth  of  the  oxide  film.  This  will  introduce 
the  presence  of  an  electric  field  in  the  oxide.  In  the  following  discussion, 

the  possible  forms  of  the  ions  will  be  explored  and  then  the  possible  elec- 

tric fields  will  be  investigated.  The  viewpoint  of  Weyl  (reference  8)  and  his 
notation  will  be  used. 

Weyl  approached  the  bonding  of  silica  from  a starting  point  of  completely 
ionized  atoms,  Si+4  and  0“2.  The  outer  electrons  on  the  oxygen  are  loosely 
bound,  and  when  brought  close  to  the  Si+4,  they  assume  more  tightly  bound 
states  which  screen  the  positive  charge  in  the  silicon.  Weyl  regarded  the 
covalent  bond  from  an  electrical  screening  viewpoint  and  indicated  this  bond 
and  the  coordination  by  writing  the  bonded  ions  adjacent  to  each  other,  i.e., 


L-2  /si\+l  nr  IVilf 
» I 0 \&  J2J  ’ or  |^4  ) ^2/ 


When  the  bonds  an  atom  is  making  are  only  partially  indicated,  only  half  or  a 
quarter  of  the  atom  is  shown. 
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If  pure  Si Oz  is  an  electrolyte,  then  Si02  must  divide  into  ions  somewhat 
like  water  divides  into  hydrogen  and  hydroxyl  ions.  There  are  two  possibili- 
ties for  the  negative  Ion;  interstitial  0“2  ions  and  the  non-bridging  oxygen 
bond.  • 


■(0;  si-  o-* 


There  are  three  possibilities  for  the  positive  ion:  an  oxygen  vacancy 


Si+4  Si+4 


++ 


a dangling  silicon  bond 


or  "triply  coordinated  oxygen" 


Energetically,  the  triply  coordinated  oxygen  is  more  favorable  because  it  re- 
sults in  better  screening,  i.e.,  the  charges  are  more  diffuse  (each  Si  would 
be  left  with  only  one  third  of  an  unscreened  positive  charge).  In  the  case 
of  the  dangling  silicon  bond,  there  is  one  extra  positive  electronic  charge. 
The  oxygen  vacancy  requires  two  positive  charges  side  by  side  and  its  exis- 
tence obviously  requires  too  much  energy. 

Diffusion  or  motion  of  the  non-bridging  oxygen  is  achieved  by  breaking  a 
bond  of  an  adjacent  oxygen  and  making  a bond  with  the  previously  unbonded 
oxygen.  Diffusion  of  the  triply  coordinated  oxygen  requires  that  one  of  the 
silicons  break  its  triply  coordinated  bond,  briefly  forming  a dangling  bond, 


Si+4 


+ 
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The  silicon  then  passes  through  the  center  of  the  triangle  made  by  the  three 
oxygpns  and  makes  a triply  coordinated  bond  on  the  other  side  by  bonding  to  a 
formerly  doubly  coordinated  oxygen.  In  the  case  of  excess  oxygen  ion  diffu- 
sion, the  extra  energy  required  to  form  the  ion  can  be  transferred  along  with 
the  non-bridging  character.  Diffusion  of  triply  coordinated  oxygen  requires 
enough  energy  to  form  a dangling  silicon  bond.  Thus,  we  would  expect  from  the 
above  that  oxygen  ions  would  diffuse  more  readily  than  silicon,  which  is  the 
observed  case  (reference  10). 

Oxide  films  can  be  grown  in  three  different  ways:  (1)  thermal,  in  an 

oxidizing  atmosphere;  (2)  anodic,  anodize  in  special  liquid  electrolytes 
(e.g.,  tetrahydrofurfuryl  alcohol  solutions);  (3)  plasma,  anodize  using  an 
oxygen  plasma  as  the  electrolyte.  During  anodic  growth,  silicon  moves  to  the 
electrolyte  Si02  interface.  It  is  also  known  that  either  light  or  another 
(doping)  source  of  holes  is  required  at  the  S i - S i 0 2 interface.  Thus,  anodic 
growth  consists  of  a half  cell  reaction 


4p  + Si  + 2Si 02  - 4 


m:  h 


and  the  transport  of  the  silicon  excess  centers  to  the  Si02-electrolyte  inter- 
face, e.g.,  in  a concentrated  nitric  acid  solution  where  oxidation  is  com- 
pleted by 

4 \(jr)3  °"2]  + 4HN03  "■  3 si °2  + 2N20 5 + 4H+ 


The  H+  ions  then  proceed  to  the  other  electrode  in  the  nitric  acid  to  com- 
plete the  circuit. 

In  the  case  of  plasma  anodization,  activation  energy  to  create  oxygen 
excess  enters  at  the  Si 02  plasma  interface  is  available  from  tne  bombarding 
oxygen  ions. 


4S i 02  ^ 02 


3e“  - 4 
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These  centers  are  transported  toward  the  Si -SI 02  interface  by  diffusion  and  the 
electric  field.  The  excess  silicon  ions  are  transported  by  the  electric  field 
toward  the  Si02  plasma  interface.  The  two  types  of  centers  can  combine  in  the 
bulk  of  the  S102. 

A thermal  oxide  could  proceed  by  a process  similar  to  plasma.  However, 
despite  Jorgensen's  experiment,  thermal  oxidation  must  proceed  mostly  by  un- 
charged oxygen  excess  centers  for  the  following  reasons. 

(a)  There  is  no  external  circuit  to  carry  the  electrons  from  the  silicon 
interface  to  the  oxygen  interface. 

(b)  Grove  (reference  10)  Indicates  that  there  are  5 x 1016  02/cc.  If 
these  were  all  uniformly  distributed  charged  oxygen  excess  centers,  i.e., 

(1  x 1017)  charges,  then  integrating  Poisson's  equation,  the  charge  density  in 

O 

a 1000  A thick  film  would  result  in  a potential  of  4.5  volts. 

(c)  It  is  known  that  in  thermal  oxidation,  the  oxygen  diffuses  toward  the 
silicon  interface  rather  than  the  silicon  toward  the  oxygen  interface.  In 
thermal  oxides,  the  negative  electrode  is  not  present,  however,  the  positive 
electrode  is  present.  Accordingly,  positive  silicon  excess  centers  can  be  cre- 
ated and  drawn  into  the  oxide  If  the  necessary  field  is  present.  Once  the  ex- 
cess center  is  in  the  oxide,  a hole  can  quantum  mechanically  tunnel  back  to  the 
silicon,  leaving  an  uncharged  center  behind.  Thus,  thermal  oxide  could  proceed 

O 

by  the  process  of  introducing  excess  silicon  centers  into  the  first  100  A of 
the  silicon  dioxide,  where  excess  oxygen  centers  combine  with  them  to  form 
stoichiometric  S102.  The  rate  of  excess  silicon  center  Introduction  would  de- 
pend on  the  field. 

In  anodic  and  plasma  oxidation,  there  is  an  applied  electric  field  present. 
According  to  Jorgensen's  experiment  (references  9 and  11),  there  is  a chemical 
EMF  present  during  the  growth  of  thermal  oxides.  In  table  1,  these  electro- 
chemical potentials  for  various  reactions  are  recorded.  The  possibility  that 
the  electric  field  present  during  oxidation  could  have  an  important  influence 
on  the  resulting  oxide  is  discussed  in  the  following  paragraphs. 

Consider  the  formula  that  Cabrera  and  Mott  (reference  12)  give  for  the 
rate  of  oxidation  of  metals 

dx/dt  = NV  v exp  [-(W+U)/kT]  sinh  [qaE/kT]  (5) 
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where:  V is  the  volume  per  metal  atom  in  the  oxide 

=*1/(2. 3 x 10+22)cm3, 

n is  the  number  of  sites  per  cm2  on  the  interface  from  which  a metal 
atom  can  be  extracted  (figure  10), 


SILICON 

°oo  O 

8888 

8888 
oooo 
Qgoo 
oooo 
o oo  o 
oooo 
2200 
oooo 

22° 
ooo 
ooo 
ooo 

888 

§88 

22° 
coo 
ooo 


OXIDE 


7 


THE  OXIDIZING  EDGE 
OF  A CRYSTAL  PLANE 


Figure  10.  The  Oxidizing  Edge  of  the  Crystal  Plane. 

v Is  the  vibrational  frequency  of  the  metal  atom  ~1012  - 1013  sec"1, 

W and  U are  the  activation  energies  of  ion  creation  and  diffusion  of 
the  metal  ion  as  shown  in  figure  11, 


THE  ION  POTENTIAL 

Figure  11.  The  Ion  Potential. 
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q is  the  electronic  charge, 
kT  is  the  thermal  energy, 

E is  the  electric  field,  and 
a is  the  distance  shown  in  figure  11. 

The  formula  can  be  applied  to  Si02  growth  data.  First,  from  Wegener  (ref- 
erence 13),  8A  of  thermal  oxide  takes  5 min  at  850°C  and  15  min  at  800°C  (also 

O O 

at  850°C,  15  min  gives  23  A and  30  min  gives  43  A).  Second,  Schmidt  (reference 

o 

14)  reports  room  temperature  growth  of  anodic  oxide  with  a growth  rate  of  25  A/ 

o o 

min  at  a constant  7 mA/cm2  and  a constant  field  of  0.26  V/A  or  3.8  A/V.  How- 
ever, there  is  a 30  V jump  during  the  first  minute  which  corresponds  to  an  un- 
known thickness,  indicating  a field  in  the  first  few  angstroms  greater  than 
the  0.26  V/A  in  the  rest  of  the  film. 

To  evaluate  this  formula,  the  following  values  were  assigned:  N = Q$s  = 

1010  cm-2;  V = 1/(2. 3 x 1022)cm3;  v = 1013  sec'1;  implying  NVv  = 2.6  x 1010 

O 

A/min.  For  Wegener's  data,  E was  obtained  by  dividing  the  value  from  table  1 

o o 

by  6 A and  solving  for  W + U,  which  was  2.36  eV  (6  A was  chosen  because  the 

o o 

rate  was  nearly  linear  from  8 A,  thus  some  thickness  less  than  8 A was  needed). 
For  Schmidt's  data,  the  above  W + U ( i . e . , 2.36  eV)  was  used  to  solve  for  E, 

O 0 0 

which  was  3.4  V/A  (30  volts  across  8.8  A gives  3.4  V/A).  Consistent  results 
are  then  obtained  with  the  above  formula  if  the  available  voltage  is  expressed 

O 

across  the  first  6 to  9 A. 

The  above  analysis  is  based  on  the  idea  that  there  is  a high  field  in  the 
first  6 to  9 A.  This  field  (0.23  to  3.4  V/A)  has  not  been  detected  after  oxi- 
dation. Since  it  corresponds  to  1013  to  1015  negative  charges  per  cm2,  it 
would  have  been  seen  if  it  was  present.  Consequently,  either  it  anneals  when 
the  wafer  is  removed  from  the  oxidizing  environment,  or  it  is  not  there.  In 
the  case  of  thermal  oxide,  an  alternative  field  is  possible.  If  each  Qss  site 
is  assumed  to  have  a positive  charge,  and  also,  opposite  it  in  the  oxide, 
there  is  an  opposing  negative  charge,  then  the  field  to  pull  the  silicon  into 
the  oxide  Is  only  present  at  the  Qss  site  and  not  everywhere,  and  the  formula 
in  equation  (5)  would  still  apply.  This  nonuniform  field  would  still  be  af- 
fected by  externally  applied  fields  so  the  above  possibilities  can  still  be 
explored  experimentally. 
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In  the  discussion,  the  form  that  ions  might  take  in  Si02  has  been  explored, 
as  well  as  the  electric  fields  that  might  be  present.  It  was  concluded  that 
the  nonbridging  oxygen  and  the  triply  coordinated  oxygen  were  the  most  likely 
form  of  the  ions  and  that  a field  produced  by  separate  individual  dipoles  is 
the  most  likely  in  thermal  oxides. 

4.  Final  Model  of  Oxide  Growth 

The  idea  of  a nonuniform  charge  sheet  leads  to  modification  of  the  Cabrera 
and  Mott  model.  Given  Qss  with  a positive  charge,  then  a nearby  negative 
charge  could  pull  this  ion  away  from  the  interface.  Where  do  the  negative 
charges  come  from?  Williams  (reference  15)  has  observed  at  room  temperature 
electron  trapping  centers  2 eV  or  more  below  the  Si02  conduction  band. 

Williams  (reference  15)  also  gives  4.25  eV  as  the  energy  difference  between 
the  silicon  valence  band  and  the  silicon  dioxide  conduction  band.  The  Fermi 
level  in  the  silicon  must  be  at  the  intrinsic  point  at  thermal  oxide  growing 
temperatures  (i.e.,  Ep  - Eva]ence  = 0.55  eV).  Then  the  above  trapping  cen- 
ters are  ~1.7  eV  above  the  Fermi  level  as  shown  in  figure  1.  Thus,  given 
these  traps  are  within  quantum  mechanical  tunneling  distance,  they  will  be 
filled  by  electrons  from  the  silicon  with  an  activation  energy  of  ~1.7  eV. 

If  these  traps  are  assumed  to  be  oxygen  excess  centers,  a complete  model  of 
the  oxidation  process  Is  obtained. 

Thermal  Oxides 

Oxygen  Is  absorbed  at  the  oxygen  Si 02  interface  forming  neutral  oxygen 
excess  centers.  The  neutral  oxygen  excess  centers  diffuse  to  the  silicon 

Interface  where  they  are  charged  negatively.  Then  they  are  drawn  toward  a Qss 

center  until  the  force  between  the  oxygen  and  the  positive  Qss  center  is  large 
enough  to  pull  a silicon  atom  out  of  the  silicon. 

Examination  of  figures  8 and  9 will  show  that  this  involves  breaking  two 

SI-SI  bonds.  The  excess  oxygen  center  can  satisfy  only  one  of  these  bonds. 

The  second  bond  becomes  an  excess  silicon  center.  The  excess  silicon  center 
diffuses  into  the  oxide  until  It  meets  and  combines  with  another  oxygen  ex- 
cess center  to  make  stoichiometric  Si02.  The  oxidation  rate  is  given  by 

£ = QssRA0Vve-(Eo-EF)/kT  (6) 
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where  AQ  is  the  oxygen  excess  center  activity 

= the  excess  center  concentration  given  as  1017cm“3  by  Grove  (reference 
10}  for  dry  oxygen  10i9  for  steam 

R is  a reaction  rate  constant 
-E  /kT 

= RQe  R'  where  ER  is  the  activation  energy  for  the  Qss  negative 
oxygen  center  reaction 

Eq  - Ep  is  the  energy  difference  required  to  charge  the  oxygen  center 
* 1.7  eV 

Grove  (reference  10)  has  shown  that  the  activation  energy  for  this  reaction 
is  2 eV;  as  shown  earlier,  the  decreasing  of  Qss  with  temperature  adds  a re- 
quirements for  an  additional  0.3  ev;  then  E^  must  be  about  0.55  ey,  half  the 
silicon  band  gap.  It  is  interesting  to  note  that  1017  centers/cm2  is  one 

o o 

center  per  210  A cube.  Thus  for  films  less  than  210  A the  basic  solubility 
must  change,  probably  increase,  causing  an  increase  in  A0and  causing  the 
initial  rapid  growth  described  in  Grove  (reference  10). 

Plasma  and  Platinum  Electrode  Oxide 

The  growth  of  plasma  oxides  and  the  platinum  electrode  oxides  of  Jorgensen 
(references  9 and  11)  differs  from  the  growth  of  thermal  oxides  in  that  the 
oxygen  excess  centers  are  charged  at  the  plasma  or  platinum  interface  and 
diffuse  or  can  be  driven  through  the  oxide  film  as  ions.  At  the  silicon 
interface,  the  negatively  charged  excess  oxygen  ion  draws  a silicon  from  the 
Qss  complex,  the  same  as  in  thermal  oxides.  The  combination  produces  a sili- 
con excess  center  which  can  be  easily  charged  positive  by  quantum  mechanical 
tunnelling  of  a hole  from  the  silicon.  The  silicon  excess  ion  is  then  drawn 
toward  the  plasma  or  platinum  Interface.  As  the  silicon  excess  ions  move, 
they  help  to  electrically  neutralize  the  space  charge  of  the  oxygen  excess 
ions  and  may  frequently  combine  with  oxygen  ions  to  form  stoichiometric  Si02. 
The  creation  of  new  oxide  in  the  middle  of  the  film  should  create  a film 
which  is  both  denser  and  more  highly  stressed  (easily  etched)  than  thermal 
oxides. 

Anodic  Oxides 

In  the  case  of  anodic  oxides,  the  silicon  must  be  drawn  from  the  Q$s  com- 
plexes by  the  applied  electric  field.  Here  the  above  formula  from  Cabrena 
and  Mott  applies.  The  positive  silicon  excess  ions  proceed  rapidly  to  the 
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electrolyte  Interface  producing  negligible  space  charge  to  modify  the  applied 
field.  At  the  electrolyte  interface,  the  silicon  excess  centers  are  oxidized 
to  stoichiometric  Si02. 

In  each  of  the  above  methods  of  oxide  film  growth,  silicon  centers  are 
left  behind  when  the  growth  process  stops.  We  would  expect  the  plasma/plati- 
num oxides  to  leave  behind  a generally  ionized  and  stressed  film.  Anodic 
oxides  might  be  very  excellent  unless  some  neutral  excess  silicon  centers  are 
created  and  left  near  the  silicon  interface. 

In  thermal  oxides,  a residue  of  dissolved  oxygen  and  excess  silicon  is 
left.  The  excess  silicon  is  next  to  the  silicon  interface  penetrating  into 
the  oxide  with  a decreasing  exponential,  the  solution  for  lossy  diffusion. 


D 


Si 


<av>  A„  Ac.  = 0 
0 Si 


A^  = Be-aX  + ceaX 


(7) 


where  a = ^ <av>  Aq/D<^ 

is  the  diffusion  constant  for  excess  silicon  centers 

A<.^  Is  the  excess  silicon  center  activity  (concentration) 

Dc.  contains  an  activation  energy  En 
si  u$i 

<ov>  Is  the  capture  constant  of  oxygen  excess  centers  for  silicon  ex- 
cess centers 

An  is  the  oxygen  excess  center  concentration  a 1/N  „ 
o ss 

B and  c are  constants  to  be  determined  by  boundary  conditions 
The  boundary  conditions  are 

1.  At  x = o the  flux  of  centers  as  determined  by  Fick's  equation  must  be 
equal  to  twice  the  rate  that  Si  atoms  enter  the  oxide,  i.e., 


JS1 


(-°S1 


VA 


Si 


\ 1 dx 

iJx=o  V eft 


Si 
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.ax 


- ce 


+aXi 


(8) 
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2.  At  x = t,  the  thickness  of  the  oxide,  the  presence  of  an  oxygen  atmos- 
phere will  oxidize  any  remaining  excess  SI  centers  thus 


Be'aX  + ce+aX  = o 


(9) 


With  the  above  boundary  conditions,  the  solution  for  ASi  is 

ASi  = K(e"aX  " ea(x'2t))/(l  + e'2at) 


(10) 


where  K = Jsi/Dsia 

Next  consider  the  worst  case  radiation  susceptibility  problem,  i.e.,  all  the 

O 

silicon  excess  centers  have  been  charged  positively  except  the  first  25  A,  (x0), 
from  which  the  charge  can  escape  by  QM  tunnelling  (reference  16).  The  amount  of 
charge  Induced  at  the  silicon  interface  and  at  the  silicon  dioxide  metal  inter- 

t . x 

face  is  equal  to  the  charge  in  the  trap  and  is  divided  — z — to  the  Si  inter- 

X ^ 

face  and  ^ to  the  metal  interface.  Thus,  the  charge  induced  in  the  silicon  is 


dx 


K e~aXo 
a (1+e-2at) 


Using  equation  (6) 


Hvl)  (£*)  - 


0 R ve-(2.3  eV)/kT 
1C  m gssVe 

a <oV> 


(ID 


(12) 


The  value  of  aQ/k  is  plotted  in  figure  12  for  x„  = 25  A and  1/a  = 12.5,  25,  50, 

a 0 

100,  200,  400,  800,  and  1600  A.  It  is  known  that  the  positive  charge  after 

0 0 o 

radiation  Is  within  200  A so  a must  be  50  to  200  A.  Notice  that  at  t = 1000  A 

O 

a variation  of  a between  50  and  400  A makes  little  difference.  Thus,  the  vari- 
ation of  hardness  with  oxygen  partial  pressure  must  be  looked  at  at  a differ- 
ent thickness. 
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5.  Interface  States 

Up  to  this  point,  no  extension  has  incorporated  Nr.  Radiation  produces 
interface  states  with  energies  distributed  across  the  forbidden  gap  of  silicon. 
Of  all  these  states,  the  most  harmful  to  the  MOSFET  function  are  those  near 
the  conduction  and  valence  bands.  Their  behavior  is  identical  to  those  which 
produce  N$s.  There  is  no  information  available  to  indicate  that  they  are 
caused  by  defects  different  than  those  which  cause  Nss.  Thus,  there  is  no 
real  extension  to  the  model  to  incorporate  Np  other  than  those  already  made 
for  Nss. 

The  production  of  by  radiation  could  be  by  removal  of  the  passivating 
oxygen  atoms  from  Si -0-Si  surface  bonds  (figure  3).  It  has  been  suggested, 
since  N is  so  easily  reduced  by  HgO,  that  there  is  some  passivation  of 
single  dangling  Si  bonds  by  OH"  ions  producing  a Si -OH  bond.  The  breaking  of 
this  OH"  bond  could  be  the  explanation  for  the  production  of  states  in  the 
center  of  the  forbidden  gap. 

The  action  which  breaks  these  bonds  could  be  direct  ionization  or,  in- 
directly, the  positively  charged  Q (particularly  the  temporary  (3  in  the 

O * l 

first  25  A),  electro-statistically  attracting  the  negatively  charged  passi- 
vating species  away  from  the  interface. 
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SECTION  IV 

SUMMARY. AND  RECOMMENDATIONS 

In  the  preceeding  discussions,  the  model  has  been  changed  and  extended  in 
the  following  ways: 

(a)  Qss  is  caused  by  an  oxygen  atom  which  is  bonded  to  three  interface 
silicon  atoms.  This  produces  an  interface  donor  state  whose  energy  is  above 
the  silicon  conduction  band. 

(b)  N$s  is  primarily  caused  by  a Si -0-Si  bond  in  which  the  oxygen  has  been 
removed  leaving  a stretched  Si -Si  bond.  This  produces  two  interface  states,  a 
donor  level  near  the  valence  band  and  an  acceptor  level  near  the  conduction 
band. 

(c)  Nr  near  the  band  edges  has  the  same  cause  as  Ns$.  The  oxygen  is  re- 
moved from  the  Si-O-Si  bond  by  direct  action  of  the  radiation  and  by  the  at- 
tractive action  of  the  positively  charged  ions  which  produce  Q . 

(d)  Qr  is  the  result  of  the  holes  which  are  produced  by  the  ionizing  radi- 
ation, being  trapped  on  silicon  excess  centers. 

(e)  Silicon  excess  centers  are  centers  where  one  oxygen  is  bonded  to 
three  silicon  atoms  in  the  bulk  of  the  Si02.  The  excess  centers  are  distrib- 
uted according  to  the  form  A^  = K(e~aX  - ea^x"^t^)/(l  + e~^at). 

(f)  The  susceptibility  to  Qr  is  proportional  to  Qss  measured  after  oxida- 
tion, before  annealing. 

(g)  Internal  electric  fields  are  important  in  oxide  growth  processes,  and 
the  growth  can  be  Influenced  by  anything  which  alters  the  field. 

The  above  are  not  an  established  part  of  the  model  but  working  hypotheses. 
It  is  unnecessary  to  perform  experiments  to  establish  each  of  the  extensions. 
Better,  use  these  hypotheses  to  analyze  current  processing  practice  and  to  pre- 
dict methods  which  would  produce  hardened  oxides.  The  test  of  these  methods 
gives  both  experimental  confirmation  and  a hardening  process  modification. 

Experiment  #1 

Aubuchon  (reference  17)  reported  that  electron  beam  evaporation  of  alumi- 
num produces  softer  oxides  than  crucible  aluminum.  This  effect  was  blamed  on 
ionization  of  the  oxide  by  x-rays  and  stray  electrons.  Others  have  reported 
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that  this  is  not  true.  Moreover,  the  effect  is  not  annealed  by  the  500°C 
aluminum  sinter,  whereas  other  researchers  have  observed  that  the  irradiate 
and  anneal  process,  if  anything,  improves  hardness  (references  1 and  4).  How- 
ever, if  there  was  irradiation  in  the  electron  beam  process  and  that  radiation 
did  charge  the  excess  silicon  centers  positive,  then  the  positive  ions  could 

o o o 

be  moved  deeper  (e.g.,  from  the  0 to  25  A region  to  the  25  A to  200A  region) 
in  the  oxide  by  an  anodizing  field  (e.g.,  applied  by  collecting  stray  elec- 
trons). Thus,  depending  on  the  particular  apparatus,  there  may  or  may  not  be 
an  effect,  which  would  explain  the  controversy.  The  experiment  then  involves 
putting  contacts  on  the  wafer  during  electron  beam  aluminum  ionization  and 
applying  both  anodic  and  cathodic  fields.  If  the  cathodic  potential  produces 

o o 

harder  devices  (by  pushing  excess  silicon  centers  from  the  25  A to  200  A re- 

O 

gion  into  the  0 to  25  A region),  then  the  ionic  nature  of  the  centers  is  con- 
firmed and  a new  hard  process  is  available. 

Experiment  #2 

An  extension  of  the  above  is  to  apply  anodic  or  cathodic  fields  during 
high  level  irradiation  (106  rads/mi n)  at  different  temperatures  to  control 
diffusion  rates  and  the  rate  of  thermal  release  of  trapped  charge.  This  could 
result  in  a modification  of  the  sintering  process. 

Experiment  #3 

The  concentration  of  oxygen  in  Si02  during  oxidation  is  5 x 1016  atoms/cc 
whereas  water  reaches  a level  of  3 x 1019  (reference  10).  Qs$  and  N_s  are 
both  much  smaller  in  steam  oxides  than  dry  oxides.  Thus,  steam  oxides  should 
be  much  harder,  which  is  not  observed.  The  answer  could  lie  with  the  expoen- 
tial  in  equations  (6)  and  (12).  If  the  activation  energy  for  charging  the 
excess  oxygen  center  is  changed,  Q is  changed.  This  could  be  caused  by  a 
variation  in  the  field.  Water  does  have  a different  vFt  (table  1)  which  means 
a different  field.  Goetzberger  (reference  18)  reports  that  fields  can  be 
applied  through  steam  during  oxide  growth  but  not  dry  oxygen,  and  such  fields 
alter  Qss.  The  experiment  is  then  to  repeat  Goetzberger's  experiment  evalu- 
ating radiation  susceptibi 1 ity.  If  the  results  prove  important,  a laser  or 
microwave  energizer  can  be  added  to  ionize  the  oxidizing  gas,  thus  permitting 
the  application  of  a field  through  dry  oxygen. 
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Experiment  #4 

The  previous  experiments  are  based  on  changing  the  location  of  the  excess 
silicon  centers  by  ionizing  the  centers  and  driving  them  about  with  electro- 
static fields.  This  cannot  be  done  with  impunity,  for  experiments  with  anodic 
growth  (references  13  and  19)  have  shown  that  such  fields  can  induce  further 
growth  and  the  introduction  of  additional  centers.  Thus  the  final  experiment 

is  to  further  explore  the  properties  of  anodic  growth.  First,  at  what  fields 

0 

does  growth  begin?  One,  perform  anodic  growth  of  very  thin  0 to  25  A films, 
continuously  measuring  thickness  versus  time  and  voltage  to  explore  the  ob- 
served 30  V initial  jump  in  applied  voltage.  Two,  induce  anionic  growth  in 
hard  oxide  films  as  a function  of  temperature.  These  two  experiments  deter- 
mine if  there  is  a charge  sheet  in  anionic  growth.  Second,  how  many  centers 
does  anodic  growth  introduce?  Irradiate  the  above  films  to  measure  Q . 

Experiment  #5 

This  is  both  an  experiment  and  reanalysis  of  old  experiments.  The  task  is 
to  gather  all  susceptibility  data  available  and  compare  equations  (6)  and  (9) 
to  look  for  a correlation. 

Experiment  #6 

O 

Evaluate  the  susceptibility  of  400  A or  less  oxides  grown  in  different 

O 

reduced  partial  pressures  of  oxygen.  Also  evaluate  at  2000  A in  reduced  and 
increased  partial  pressure.  This  experiment  should  use  etch  back  or  some 
other  means  of  locating  the  charge  centroid. 

The  above  experiments  and  tasks  are  then  recommended  to  verify  or  reject 
the  model  extensions  and  to  investigate  possibilities  for  hardened  processes 
and  hardness  assurance  screens. 
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